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Abstract 

A non-Abelian version of the Hirota-Miwa equation is considered. In an earlier paper [Nimmo 
(2006) J. Phys. A: Math. Gen. 39, 5053-5065] it was shown how solutions expressed as quasideter- 
minants could be constructed for this system by means of Darboux transformations. In this paper 
we discuss these solutions from a different perspective and show that the solutions are quasi-Pliicker 
coordinates and that the non-Abelian Hirota-Miwa equation may be written as a quasi-Pliicker rela- 
tion. The special case of the matrix Hirota-Miwa equation is also considered using a more traditional, 
bilinear approach and the techniques are compared. 

1 Introduction 

There has been interest recently in noncommutative version of some of the well-known soliton equations 
including the KP and KdV equations [1-10]. In many cases, the noncommutative version of the equation 
is obtained by considering the compatibiUty of the same Lax pair that is used in the commutative case. 
The only difference between the commutative and noncommutative cases being that the assumption that 
coefficients in the Lax pair commute is relaxed. 

The Hirota-Miwa equation, or discrete KP equation, [11, 12] is the fully discrete three dimensional 
integrable system 

(a2 - a3)r iT 23 + (as - ai)r 2T,3i + (ai - a2)T3Ti2 = 0, (1) 
where Ofc are constants, r — T(ni,n2,"-3) and subscripts j denote a forward shift in the corresponding 
discrete variable ni, n2 or n^. This is the compatibility conditions for the linear system 

- '/>,j = U,j(t>, (2) 

for i < j = 1, 2, 3 where 

Uij = (flj - Qj)^^. (3) 

If one considers the linear system ([5]) in a noncommutative setting, the parameterisation ^ is not 
appropriate but nonetheless we obtain a nonlinear system for the Uij which we call the non-Abelian 
Hirota-Miwa equation. This system was derived in [13] and solutions in terms of quasideterminants [14] 
were derived by means of Darboux transformations. 

It has long been known [15] that equations in the KP hierarchy can be interpreted as Plucker relations 
on an infinite dimensional Grassmann manifold. In this paper we establish a result which hints at a 
corresponding situation in the noncommutative case as we show that the quasideterminant solutions are, 
in the terminology of Gelfand et al [14,16], quasi-Pliicker coordinates and the non-Abelian Hirota-Miwa 
equation is a quasi-Pliicker relation. 

This paper is arranged as follows. In Section 2, the properties of quasideterminants that are used 
in the rest of the paper are described. In Section 3, the main results from [13] are summarised and 
the quasideterminant solutions found there are presented. The next two sections give two contrasting 
approaches to direct verification of these solutions. In the first, applicable in the general non-Abelian 
case, the solutions are recognised to be quasi-Pliicker coordinates and the direct verification is effected 
using quasi-Pliicker relations and other properties. Finally, we consider an alternative approach which 
works only for the matrix Hirota-Miwa equation. In this it is shown how standard Pliicker relations are 
used to verify the solutions. 
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2 Preliminaries 



In this section we will state the properties of quasideterminants that are used in this article. The reader 
is referred to the original papers [14, 16] for a more detailed description and proofs of the results described 
here. 



2.1 Quasideterminants 

For an nxn matrix M = (aij) over a ring TZ (noncommutative, in general) there are quasideterminants 
written as \M\ij for i,j = 1, ... ,n which are also elements of TZ. They are defined recursively by 



(4) 



In the above represents the ith row of M with the jth element removed, Cj the jth column with the 

ith clement removed and A'P'^ the submatrix obtained by removing the ith row and the jth column from 
A. Quasideterminants can also denoted as shown below, by boxing the entry about which the expansion 
is made 

M''3 c:- 



|M| 



In the case n 
For n = 2, if M = 



1, let M = (a), say, and then there is just one quasideterminant |Af|i^i = | | a \ 
a b 



a. 



c d 

\M\lA 



, then there are four quasideterminants 





c 



\M\ 



2,1 



a b 

From this we can obtain the matrix inverse 



a — bd ^c, 
c — db~^a, 



\M\l,2 
\M\2,2 



d 



b — ac ^d, 
d — ca~^b. 



A^-i _ A« - M-^c)-! (c - db-'^a)-'^ 
~ {{b - ac-^d)-^ {d-ca-^b)-^ 

which is then used in the definition of the 9 qiiasidc;terniinants of a 3 x 3 matrix. 

To explain results in the simplest way, we will sometimes write M as a block matrix 



M 



A B 
C d 



where d € TZ, A is a square matrix over TZ of arbitrary size and B, C are column and row vectors over 
TZ of compatible lengths. For this matrix, 

A B 



C d 



d-CA-^B. 



(5) 



2.2 Commutative cases 

If the entries in M commute then the above becomes the familiar formula for the inverse of a 2 x 2 
matrix with entries expressed as ratios of determinants. Indeed this is true in general; if the entries in 
M commute then 



r ■ 



\M\ 



(6) 



the alternating sign arising from the reordering of rows and columns in the numerator. 

In the special case that the ring TZ is the (noncommutative) ring of N x N matrices over another 
commutative ring, the (z, j)-quasideterminant of M, \M\ij is itself a,n. N x N matrix with (A;,^)th entry 



{\M\ij)k,i = 



k.l 



(7) 
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2.3 Invariance under row and column operations 

The quasideterminants of a matrix have invariance properties similar to those of determinants under 
elementary row and column operations applied to the matrix. Consider the following quasideterminant 
of an n X n matrix; 













EA EB 












n,n 


FA + gC FB + gd 


n.n 



: g{d-CA-^B) 



(8) 



This formula describes the effect on the quasideterminants of a matrix of performing elementary row 
operations, involving left multiplication, on the matrix. The zero block in the first matrix means that 
operations which add left-multiples of the row containing the expansion point to other rows are excluded 
from consideration, but all other elementary row operations are included. We see that left-multiplying 
the row containing the expansion point by g has the effect of left-multiplying the quasideterminant by 
g and that all other operations leave the quasideterminant unchanged. There is analogous invariance 
under column operations involving multiplication on the right. 



2.4 Noncommutative Jacobi Identity 

There is a quasideterminant version of Jacobi's identity for determinants, called the noncommutative 
Sylvester's Theorem by Gelfand and Retakh [16]. The simplest version of this identity is given by 

ABC 

^/3 = z.m-z.mnnn nm- (9) 

E h \i\ 



2.5 Quasi-Pliicker coordinates 

Given an (n -1- fc) x n matrix A, denote the ith row of A by Ai, the submatrix of A having rows with 
indices in a subset / of {1, 2, . . . , n + fc} by A/ and ^{i,...,n+fc}\{i} by A^. Given i,j £ {1, 2, . . . , n -t- fc} 
and / such that #/ = n ~ 1 and j ^ /, one defines the (right) quasi-Pliicker coordinates 



Ai 




Ai 


A, 


71 S 





Ai 

A, 
A, 







1 



(10) 



for any column index s £ {1, . . . ,n}. The final equality in (fT(I|) comes from an identity of the form Q 
and proves that the definition is independent of the choice of s. 
The following are easy consequence of the definition: 



if i e /, (and is not defined if j were in /) 
1 



ij jk ik' 

Remark. A useful consequence of (fTQ|) and (fT3]l is the identity 



A' 





~i 


A' 





A' 







A' 


1 


A> 


1 




Ai 






which shows that quasideterminants of this form may be inverted very simply. 

There are also the following properties which are proved in [16]: skew- symmetry 

ij jk ik 1 

for any {«, j, fc} £ / and quasi-Pliicker relations 

for any i ^ M, where L, M C {1, . . . , n -I- fc} such that #L - 1 = #Af = n — 1. 



(11) 
(12) 

(13) 
(14) 



(15) 



(16) 



(17) 
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3 A non-Abelian Hirota-Miwa equation 

This section is mainly a review of the results found in [13]. When discussing lattice equations one 
generally thinks of the dependent variable as a real or complex valued function of one or more integer- 
valued variables. Here we wish to abstract this notion as far as possible and think of a lattice in a more 
general way which encompasses the usual types of lattice systems. Similar ideas were also explored by 
Matveev in [17]. We suppose that dependent variables take values in a ring TZ and the lattice is defined 
by means of mappings Si : TZ —>■ TZ, i € N, which we think of as moving between neighbouring sites on 
the (possibly infinite dimensional) lattice. The particular examples we have in mind for these mappings 
include the usual shift operators Si{(j){ni)) = (j){ni + 1) or q-shift operators Si{(j){ni)) — (j){qini), but 
also Darboux or Backlund transformations Si{(f>) — (f> mapping between different solutions of the lattice 
equations. 

With these sorts of examples in mind, we require that the Si have the following properties: 

1. all Si are invertible and mutually commute, so that motion on the lattice is reversible and path 
independent 

2. each Si is a linear homomorphism so that the same (linear and nonlinear) relations satisfied at one 
lattice site are satisfied at all lattice sites, 

and we refer to such mappings as lattice operators. We will also use the shorthand notation Xj = Si{X) 
and X, = 5-1 (X). 

Now consider lattice operators s^, (i G N), and invertible Uij G TZ and the system of linear equations 
for (f> 

0,^-(/)j+C/y(/):=s,((/))-Sj(?!)) + C/y(/) = O, (i,ieN). (18) 
This system is algebraically over-determined and so for any k, 

Urj + Ujk + Uk^ = 0. (19) 

This includes degenerate cases which arise when some of i, j, A: coincide, namely Uu = and Uij + Uji = 0. 
The compatibility conditions (j)^ij — <j)ji give the nonlinear conditions 

Uij^kUik — UikjUij (20) 

for i,j,k distinct (if any of i,j,k coincide this is a empty condition) together with 

Uij^k + UjkA + Ukij = 0. (21) 

By introducing Vij = Uij^j-^, one may rewrite (j2ip as the "dual" equation 

V,j + Vjk + Vk^ = 0, (22) 

and one may also rewrite (j20p in terms of Vij as 

V^,Vk,j=VkV,^-^. (23) 

It is straightforward to show that the Uij satisfy (I19p and (|20|) if and only if the Vij satisfy ([22|) and 
dU) and so it suffices to consider either the equations for Uij 

Uij + Ujk + Um ^ 0, U,j,kU,k = U^k,jU,j, (24) 

or for Vij 

Vj + Vjk + Vm - 0, VkV,j-k = VjVkj. (25) 

We call either of these equivalent systems the nonahelian Hirota-Miwa equation. 

If we fix on one particular three dimensional sublattice (defined by si, S2, S3 say), the linear equations 

J +C^u<^ = 0, (z,j e {1,2,3}) (26) 

are compatible provided ([M]) or (P5|) are satisfied for A;} C {1,2,3}. The remaining dimensions in 
the lattice defined by S4, S5, • • • = di, d2, . ■ . can be interpreted as Darboux transformations. In [13] these 
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Darboux transformations were exploited in order to construct solutions of the non-Abelian Hirota-Miwa 
system expressed in terms of quasideterminants. Finally in the section, we will state these results. 
Making the ansatz 

Uij = XijUijX^ ^, (27) 

where Xi G TZ are invertible and G TZ, the nonlinear conditions (j20p are satisfied for any Xi provided 
the Qij themselves satisfy this condition. We also take an = so that Uu = 0. With these choices, ([25)1 
are solved by writing 

V,, = Yr^a,,Y,j, (28) 
where Yj, — X^^ . In the commutative case, the skew-symmetry conditions 



— Uij + Uji — XijaijXj^ ^ + Xj^iGjiX- ^ 



X, 



X, 



' X, 



can be solved identically by the additional ansatz Xi = r^i/r and taking + aji = 0. The remaining 
equation in (12411 now becomes the standard Hirota-Miwa equation ([T]) and it is for this reason that we 
refer to (|24p and as non-Abelian Hirota-Miwa equations. 

Now let Xi, i — 1,2,3, and Uij i,j e {1,2,3} construct, through (P7| . a vacuum solution Uij of 



(PSj) . For example, one could choose Xi = X2 = X3 = 1 and 
eigenf unctions, 



, 9n i.e. solutions of 
one may construct the row vector 9 ~ {9i 



ai — OLj. Now choose n vacuum 
with the vacuum choice for Uu made above. From these 



. , 6'„) satisfying 



s,{9)-s,{9)^U,,9^Q. 



(29) 



Then, after n applications of Darboux transformations, the solution obtained is Uij — XijUijX^ ^ 
where [13] 



X, 



S^{9) 



X,.. 



(30) 



4 Direct verification of solutions using quasideterminants 

In this section we use the properties of quasideterminants stated in ^Uto effect a direct verification of the 
solutions described in pop . By recognising that the quasideterminant in this equation is a quasi-Pliicker 
coordinate pop . it may immediately be inverted using (jlSp to obtain 



x; 



X- 



S^{9) 



(31) 



Then, use the linear equation and the invariance properties ([5]), we get 



S^i9) 



-\9) 
-\9) 



-,n — 3 



s7 -{9) 



(9) 



(32) 



In fact, this is the first key step in the direct verification, we can go one step further with Vij. By 
making use of (j29p and ([SJ we may write all but three rows of the quasideterminants in a way which is 



5 



independent of the particular choices of i and j. Let X^™) denote s™(X) for any choice of fc = 1, 2 or 3. 



using 


dH]) and 


® 


that 








1 











n 

(7 






















n 




n(n — 3) 













Q(n-2) 







6)(n-3) 














6)(n-2) 


1 



(33) 



To complete the direct verification of these solutions we must show that they satisfy 

V^^ - 0, V,, + Vy, = 0, + V^jfe + T4, = 0. 



(34) 



One consequence of ([8]) is that a quasideterminant in which the row containing the expansion point 
is the same as another row is zero. The second factor in Va is such a quasideterminant and so Va — 0. 
Next define the (n + 3) x n matrix 

^.2 



e 



^,3 
Q(n-2) 



Q{n-3) 



which will play the role of the matrix A in the definition of quasi-Pliicker coordinates (|T0|) . and let 
/ = {6, . . . , 77 + 3} so that 



0in-3) 



Then the quasideterminants appearing in p3p can be expressed as quasi-Pliicker coordinates, 

6* 6* 



/U{4} 
'5,i 



(6) 



6»("-2) 



1 



(9) 



(n-3) Q 



(n-2) 



Hence, 



v.. - /U{4} /U{j} 
— '5,i 'i,4 ' 



(35) 

where here and from now on we omit reference to the matrix O. The properties (|14p . p6p and p7p now 
take the particular form 



Ill{d} Ill{d} ^ IU{d} 
a,b b.c a,c 



IU{c} Ill{a} _ _ 

' a,b ' b,c ~ 

IU{c} IU{d} TU{b} IU{d} _ -I 

' a, fa ' b,a ' a,c ' c,a — 



r,/U{6} 



(36) 
(37) 

(38) 
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Using ([M ]) and (P7 | . 

V-.A-V--- /U{j} /U{4} /U{z} 

_ /U{4}. /U{i} /U{4} /U{i}x 

_ „-fU{4}. /U{i} /U{J}^ „ 
— '5,1 V i,4 ' iA ) — ^■ 

Similarly, using the quasi-Pliicker relation ([55)1 as well as p6p and ([57)1 , we have 

y,, + fl, + V^,, - /U{4} /U{2} /U{4} /U{3} /U{4} /U{1} 
fl2 -I- V23 + ^31 — ' 5,1 '1,4 +'5,2 '2,4 +'5,3 '3,4 

_ /U{4}. /U{2} /U{4} /U{4} /U{3} /U{4} /U{1}^ 
~ '5,1 Vl,4 1 '1,3 '3,2 '2,4 +'1,3 '3,4 J 

_ ^■fU{4}/ /U{2} /U{4} /U{2} /U{3}^ 

— '5,1 l'l,4 '1,3 '3,4 '1,4 ) 

_^-fU{4}.. /U{4} /U{2} /U{3} /U{2}^^ /U{2} _ „ 

— '^5,1 U ~ ^^1,3 ^^3,1 ~ ^1,4 M,l Fl,4 — 

and so the verification is complete. 

This verification gives a suggestion that quasi-Pliicker relations play a similarly important role in 
direct verification of solutions do in the noncommutative situation as standard Pliicker relations do in 
the commutative case and we intend to pursue this idea more fully in future work. 



5 An alternative method for the matrix Hirota-Miwa equation 

In the case that TZ is the ring of r x r complex matrices, an alternative, more familiar direct approach is 
available. This involves expressing the matrix entries as ratios of determinants, reexpressing the system 
()25p in bilinear form and verifying the solutions using identities arising from Laplace expansions. In 
this section, we will contrast this more familiar, bilinear approach with the direct verification using 
quasideterminants described in the last section. We will see that the proof a similar sequence of steps 
but, in contrast with the proof involving quasideterminants, has a rather ad hoc feel to it, in which the 
details of the calculation such as keeping track of the correct power of —1, are rather intricate. 

In the case that we seek solutions Vij of ([^ has solutions in this matrix ring, the ingredients in 
the solutions ([55)1 . the elements Ok are themselves also (invertible) r x r matrices and the concatenation 
of these, 9 is a, r x rn matrix. This expression for the solutions p3p as quasideterminants can then be 
reexpressed in terms of a product of matrices whose entries are ratios of rt x n determinants, 

V^, - $ ^, (39) 



where for each i, j, Fi is scalar and Gi and hij are r x r matrices. For a more compact notation, allowing 
determinants to be expressed in terms of columns rather than rows, 9 is replaced by its transpose so that 
now ^ is a rn X r matrix. The specific expressions are 

F, = \st\9) 9 ^(1) ■•• (40) 

and the (p, q)th entries of Gi and hij are 

(G,)p, = (-l)™+^+«|s-i(0)^ 9 0(1) ••• 0("-i)|, (41) 

and 

{h,,),, = {^lY+<'\s-\9), s-\9) 9 ... (42) 

respectively, where 9p denotes the p"^ column of 9 and 9q means all the columns of 9 except the g**^. 
If Vij is to satisfy p4p then we must show that 

h,, = 0, G,h,, + Gjhj, = 0, Grh^jFk + GJh,kF^ + GkhkrFj = 0. (43) 

The first of these is obvious since if j = i then two columns in ha arc the same. We will now show how 
the other two conditions may be verified. 
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For each p, q consider the 2ni x 2rn determinants 



Lij — 



9 •• 






• • 







•• 








s^\B) 6 ■■ 


gin-3) Q(n-2} 





The expansion of Lij by rn x rn minors is, up to a sign independent of i and j, 



^ ^ )pm (^ij )mq — {Gihij^^ 



m—1 

By using elementary row and column operations we also see that Lji = ~Lij, and so we have verified 
the second conditions 

Gihij + Gjhji = 0. 

For the final condition, it is necessary to introduce auxiliary matrix variables Jy where 
By a similar expansion, the determinants 



(44) 
(45) 





Q{n-2) 




• • 





• • 







9 ■■■ 



This determinant can also be rewritten as 
So we get 

GjJij — GiFj, 

for any 

Finally, consider expansion of the determinants 



sY\9)^ 9 ... 


g{n-2) 


Up 


. . 





... 





Up 


s-\9) e .. 


g{n-2) 



— ( ly {GjJij)pq. 



(46) 



s-'9 9 .. 


g{n-3) n{n-2) 


s]:'9 {sT^9)p 


(0("-2))^ 


. . 





. . 





sl^9 {sT^9)p 


(0("-2))^ 


... g{n-3) 6)("-2) 



= {"Jkjhki + hjiFk — Fihjk)pq. 
After elementary row and column operations, this determinant can also be show to vanish and so we get 

- Jkjhu + hy,Fk - hjkF^ = 0. (47) 
Using (gH) and dH]) and then (gT]) it follows that 

GihijFk + GjhjkFi + GkhkiFj = Gj{— hjiFk + hjkFi + Jkjhki) = 0. 
This completes the alternative direct verification in the matrix case. 



6 Conclusions 

In this paper we have considered two direct approaches to verifying quasideterminant solutions of a 
non-Abelian Hirota-Miwa equation found in [13]. In the first approach, the solutions are expressed in 
terms of quasi-Pliicker coordinates and the verification is achieved by using properties of these objects, 
including quasi-Pliicker relations. The second approach, available only in the matrix case and not 
involving quasideterminants, uses an expression for the solutions as matrix of ratios of determinants. In 
the first approach, the machinery used is more sophisticated but the verification is quite straightforward 
and applicable whatever the nature of the non-Abelian equation, matrix or quaternion or otherwise. The 
second verification needs less sophisticated tools (transformation to a bilinear form and standard Pliicker 
relations) and is not always available since determinants are not defined in the general non-Abelian case. 
A direct quasideterminant approach has also been applied to noncommutative versions of the KP and 
mKP equations [10, 18] and it is hoped that this can be developed into a more widely applied direct 
approach to noncommutative integrable systems. 
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